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To analyze the crucial role of fluctuation and relaxation effects for the function 
of the human brain we studied some statistical quantifiers that support the infor- 
mation characteristics of neuromagnetic brain responses (magnetoencephalogram, 
MEG). The signals to a flickering stimulus of different color combinations have been 
obtained from a group of control subjects which is then contrasted with those of a 
patient suffering photosensitive epilepsy (PSE). We found that the existence of the 
specific stratification of the phase clouds and the concomitant relaxation singular- 
ities of the corresponding nonequilibrium dynamics of the chaotic behavior of the 
signals in separate areas in a patient provide likely indicators for the zones which 
are responsible for the appearance of PSE. 

PACS numbers: 05. 45. Tp; 87. 19. La; 89. 75. -k 



I. INTRODUCTION 

The study of dynamical time series is gaining ever increasing interest and is applied and 
used in diversified fields of natural sciences, technology, physiology, medicine and economics 
[3, 0, 0, 0, 0, 0, 0, 0], to name only a few. The majority of natural systems can be consid- 
ered dynamical systems, whose evolution can be studied by time series related to relevant 
variables on a suitable time scale. These series are often characterized by a pronounced time 
and spatial synchronization or coherence, chaotic or robust behavior. 



*Electronic address: rmy@theory.kaz an~spu.ru| 



2 



When analyzing time series data with linear methods, one can follow certain standard 
procedures. Moreover, the behavior may be described by a relatively small set of parameters. 
For a nonlinear time series analysis, this is not necessarily the case. While standardized algo- 
rithms exist for the analysis of the time series data with nonlinear methods, the application 
of these algorithms requires considerable knowledge and skills on the part of the user. 

In a nonlinear time series analysis one starts out with a reconstruction of the state 
spaces from the observed data (3, @, GU 11, Although the embedding theorems 111 
provide an important means of understanding the reconstruction procedure, likely none of 
them is formally applicable in practice. The reason is that they all deal with infinite, noise 
free trajectories of a dynamical system. It is not obvious that the theorems should be 
"approximately valid" if only the requirements are "approximately fulfilled", for example, if 
the data sequence is long, although finite and not completely noise-free. 

A possible way to study the manifestation of physical properties of random processes 
(and the Markov random processes (MRP) in particular) in time series originates from 
the theory of nonequilibrium statistical physics. The history of the fundamental role of 
stochastic processes in physics dates back a century to the Markov representations [3] of 
random telegraphic signals. Such processes still find application in models of contemporary 
complex phenomena. A few typical examples of complex physical phenomena modeled by 
the Markov stochastic processes are: kinetic and relaxation processes in gases [l5j] and 
plasma [l6|, condensed matter physics (liquids [13], solids (HI, and superconductivity fl9j). 
astrophysics [2(|, nuclear physics [111, for certain quantum relaxation dynamics j22| and in 
classical [23[ physics. At present, we can make use of a variety of statistical methods for the 
analysis of the Markov and the non-Markov statistical effects in diverse physical systems. 
Typical examples of such schemes are the Zwanzig-Mori's kinetic equations [2J], generalized 
master equations and corresponding statistical quantifiers |25| , the Lee's recurrence relation 



method [26|, the generalized Langevin equation (GLE) [271 ]. etc. 

In this paper we shall study the crucial role of relaxation and kinetic singularities in brain 
function of healthy physiological and pathological systems for the case of photosensitive 
epilepsy (PSE). In particular, we seek marked differences in large space and times scales 
in the corresponding stochastic dynamics of discrete time series that could in principle 
characterize pathological (or catastrophic) violation of salutary dynamic states of the 
human brain. As a main result, we show here that the appearance of distinct differences 
in the relaxation time scales and extraordinary stratification of the phase clouds in the 
stochastic evolution of neuromagnetic responses of the human brain as recorded by MEG 
may yield evidence of pronounced zones responsible for the appearance of PSE. 



II. STRATIFICATION IN THE PHASE SPACE AND STOCHASTIC 
PROCESSES IN COMPLEX SYSTEMS 

The phase space plays a crucial role in determining the singularities of stochastic dynamics 
of the underlying system. A set of the dynamical orthogonal statistical variables describing 
the dynamical state of the complex system is a feature important in a proper construction of 
the phase space and analysis of the underlying dynamics. Let us consider an k-dimensional 
vector of the initial state A° k = (xi,X2, x%, ....Xk) and an k-dimensional dynamic vector of the 
final state A^ +m = {x m+u x m+2 , x m+3 , ....x m+k ) , where k + m = N, k,m = N, N — 1, N — 
2, N/2 — 2, N/2 — 1, A^/2 and N denotes the sample length. From the discrete equation 
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of motion 



Axi Xl+1 Xl -{A-l} Xl 



At t i+1 -ti r 

U+i - U = T, (1) 
we obtain the equation of motion of the dynamical vectors of state A™ +fc as: 

AA™ 

— t-^— = HA™, 
At 3 

i2=i(A-l), (2) 
r 

where the shift operator A acts as Axj = Xj+i and j = m+k. 

By applying successively the quasioperator L to the dynamic variables A m {t), t=mr, 
where r is a discrete time step, j = m + k, we obtain the set of dynamic functions B„(0) = 
L A°(0), n > 1. By using the variables B n (0) one can find a formal solution of the evolution 
of these equations in the form 

A^ +fc (mr) = {l + ztL\ A° fe (0) = ^ f } _ B^. (3) 
k J j =0 3\ m J)- 

However, the use of this structure is generally not the most convenient one. An advan- 
tageous approach consists in the use of the orthogonal vectors W n , as detailed below, by 
use of the Gram-Schmidt orthogonalization procedure [28[ in place of the set of variables 
B n (0). Thus, we work with this new set of dynamical orthogonal vector variables, where 
the average < ... > should be read in terms of scalar products and 5 n ^ m is the Kronecker's 
symbol. In doing so, we find the the recurrence formula in which the "senior" variables at 
order "n" W n = W n (t) are related through with the "junior" variables of order m < n; i. e. 



|29|,|3fl 



W = A°(0), Wi = {iL -Ai}W 0) ... 
(W n W m ) = 5 n/m (|W n | 2 ), 
W n = {zL-A n }W n _ 1 + A n _ 1 W„_ 2 + ..., n>l. (4) 

Here, we have used the kinetic parameters given by the Liouville's quasioperator L as follows 
I!: 



(W n LW n ) 

n (|w n p) ' 

A (|W n | 2 ) 

where A n = with the parameter Q n denoting a general relaxation frequency. The set 

of frequencies A n describes the spectrum of the Liouville's operator L. 

Besides these relaxation time measures, we in addition shall consider information mea- 
sures that are based on time correlation functions which assume the role of memory functions 



Mi(t),i = 0, 1, 2, see in Refs. [30(; i. e., 
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Note that the quantity with i 
the vector A£. 

III. ANALYSIS OF TIME SERIES FOR THE EXPERIMENTAL DATA OF PSE 

Now we can proceed directly to the analysis of the experimental data: the MEG signals 
recorded in a group of nine healthy human subjects and a patient with (PSE) [3lJ. PSE is 
a common type of stimulus-induced epilepsy, defined as recurrent convulsions precipitated 
by visual stimuli, particularly flickering light. The diagnosis of PSE involves the detection 
of paroxysmal spikes on an EEG in response to the intermittent light stimulation. To elu- 
cidate the color- dependency of PS in normal subjects, brain activities subjected to uniform 
chromatic flickers with whole-scalp MEG have been measured in Ref. [3l|). 

The same subjects and the data set were part of an earlier study in Ref. [3l[; however, 
we shall mention the relevant details for the sake of completeness. Nine healthy adults 
(two females, seven males; with age ranging from 22-27years) voluntarily participated. Two 
additional age-matched child control subjects, and one more photosensitive patient (age 
14yr) under medication (sodium valporate), were also studied. All subjects were right- 
handed and were explicitly informed that flicker stimulation might lead to epileptic seizures. 
They gave their written informed consent before recording. The subjects were instructed 
to passively observe visual stimuli with minimal eye movement. During the testing session 
for the photosensitive patient, pediatric neurologists were present to monitor their health 
condition precautionary measure. 

The subjects were screened for photosensitivity and personal or family history of epilepsy. 
The experimental procedures followed the Declaration of Helsinki and were approved by the 
National Children's Hospital in Japan. The stimuli were generated by the two video projec- 
tors and delivered to the viewing window in the shield room through an optical fiber bundle. 
Each projector continuously produced a single color stimulus. Liquid crystal shutters were 
located between the optical device and the projectors. By alternative opening one of the 
shutters for 50 ms, 10 Hz (square-wave) chromatic flicker was produced at the viewing dis- 
tance of 30 cm. Three color combination were used : red-green (R/G), blue-green (B/G), 
and red-blue (R/B). CIE coordinates were x=0. 496, y=0. 396 for red; x=0. 308, y=0. 522 
for green; and x=0. 153, y= 0. 122 for blue. All color stimuli had the luminance of 1.6 
cd/m 2 in otherwise total darkness. In a single trial, the stimulus was presented for 2s and 
followed by an inter-trial interval of 3s, during which no visual stimulus was displayed. In a 
single session, the color combination was fixed. 

Neuromagnetic responses were measured with a 122-channel whole-scalp neoromagne- 
tometer (Neuromag - 122; Neuromag Ltd. Finland). The neoromag-122 has 61 sensor 
locations, each containing two originally oriented planner gradiometers coupled with dc- 
SCUID (superconducting quantum interference device) sensors. The two sensors of each 
location measure two orthogonal tangential derivatives of the brain magnetic field compo- 
nent perpendicular to the surface of the sensor array. The planner gradiometers measure 
the strongest magnetic signals directly above local cortical currents. From 200 ms prior 
responses were analog-filtered (bandpass frequency 0.03 - 100 Hz) and digitized at 0.5 kHz. 
Eye movements and blinks were monitored by measuring an electro-oculogram. 



(W i( 0)W i( t)) 

Ivm (| Wi (o)| 2 ) • 

corresponds to the temporal autocorrelation function of 
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Figure 1: The single phase planes of the phase portrait of the patient with PSE for MEG signals 
from an R/B combination of the light stimulus : a) plane (Wo(i), Wi(i)); b) plane (Wo(t),W2{t)); 
c) plane (W (t), W 3 (t)); d) plane (W 1 (t),W 2 (t)); e) plane (W 1 (t), W 3 (t)); f) plane (W 2 (t), W 3 (t)). 



The trials with MEG amplitudes > 3000 fT/cm and/or electro-oculogram amplitudes 
> 150// V were automatically rejected from averaging. The trials were repeated until > 
80 responses were averaged for each color-combination. The averaged MEG signals were 
digitally lowpass- filtered at 40 Hz, and then the DC offset along the baseline (—100 to ms) 
was removed. At each sensor location, the magnetic waveform amplitude was calculated as 
the vector sum of the orthogonal components. The peak amplitude were normalized within 
each subject with respect to the subject's maximum amplitude. The latency range from 
— 100 to —1100 ms was divided into 100 ms bins. Then, the peak amplitudes were calculated 
by averaging all peak amplitudes within each bin. It would be important to mention that 
no clinical photosensitive seizures were induced during the experiment. This also confirms 
the better detection power of this analysis than normal seizure detection procedure. 



IV. RESULTS AND DISCUSSION 



Here we present the information-theoretic analysis for the presence of PSE, based on 
the time behavior of the dynamical variables and the stratified structure of phase spaces. 
Some results of our quantities as derived from the theory presented here, are depicted in 
Figs. 1-7. Our results for nine healthy subjects and for the patient with PSE in comparison 
are shown in Figs. 8 - 10 . Among them are: the time traces of the MEG signals (Wo), 
and for the three junior dynamical orthogonal variables (Wi), i = 1, 2 and 3; 2), the phase 
space created by the points with coordinates (Wj), i = 0, 1, 2 and 3; 3) the phase space, 
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Figure 2: The time dependence of the first four orthogonal variables for MEG signals in healthy 
in No. 3 ( sensor m = 13, for an R/ B combination of the light stimulus: for healthy subject No. 
4 : a) W (t); b) W x {t)- c) W 2 (t);d) W 3 (t) for an R/ B combination of the light stimulus. 



filled by the trajectories (Wi(t)), i = 0, 1, 2 and 3; 4) the time dependence of the first four 
quantifiers: the time correlation functions (TCF) M (t) and the first three junior memory 
functions (MF) Mj(t), for i — 1,2 and 3. The results of the experiment for a red-blue (R/B) 
and a red-green (R/G) combination of color signals are used in all of the figures. 

As an example, the analogous results for the patient with PSE(sensor No. 10) are pre- 
sented in Fig. (1). The obtained results possess the clearly visible inconsistent character. 
The comparison clearly shows a characteristic behavior of the dynamic variables (Wi(t)) (i 
— 1,2 and 3). The time dependence of the variables {W\(t)) presents the time behavior 
of the orthogonal velocity of the signal recording in discrete form. The next-order dynamic 
variable (W^t)) describes the orthogonal acceleration, the variable (Wsit)) depicts the lon- 
gitudinal orthogonal energy current, etc. The signals (Wi(i)) in the patient with PSE can be 
characterized as regular noise. The phase clouds formed by manifold phase points drastically 
differ from the healthy ones in comparison with the patient with PSE, see, Fig. 1. 

The stratification of the phase clouds and the existence of the stable pseudoorbits are 
more visible for the healthy . In the patient with PSE (see, Fig. (1)) the phase stratification 
disappears. The phase clouds can be characterized by symmetrical nuclei, they have spatial 
homogeneity. The phase trajectories for the healthy form broken lines. 

For the patient with PSE the pictures of the phase trajectories contrast sharply with the 
case of the healthy. The phase trajectories are packed tightly within the restricted areas 
of the phase space. The drastic difference in the typical scales of the dynamic variables 
(Wi(t)) and in the size of the phase space for the healthy and for the patient with PSE 
(Fig. 1) are striking. This difference ranges from 3 times (for (Wo(t)) ) to 10 times(for 
(Wi(t)),i = 1,2, and 3) and from 10 times (for the phase plane (W (t), Wi(t)) to 80 times 
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Figure 3: The phase portraits for MEG signals in healthy person No. 3 (sensor number m = 13) 
for an R/ B combination of the light stimulus. The phase space has been created by the phase 
points F itj (t) = (Wi(t), Wj{t)), 1 < i + j < 3. 

in the corresponding clouds of phase planes (W (t), Wi(t)), i = 2, 3 and (Wi(t), Wj(t)) with 
1 <M<3. 

Thus, the signals in the patient with PSE with sensor number 10 differ from the healthy 
subjects consists in the drastic change of the fluctuation scales of the dynamic orthogonal 
variables and the space sizes of the phase clouds. The similar difference of the scales con- 
stitutes the values from 10 to 80 times. This observation let us note the specific role and 
behavior of the sensor with number n=10 in the formation of PSE mechanisms! The differ- 
ence in the scales of the orthogonal dynamic variables and in the sizes of the phase clouds 
is more drastic for sensors with numbers n = 10, 5, 23, 40 and 53. 

From the time dependence of the initial TCF Mo(t) and the first three MF's of the junior 
order Mj(t), i = 1, 2 and 3 it is possible to see great difference in the behavior of the time 
functions for the healthy and for the patient with PSE. One can observe a large-scale time 
correlation in the healthy subjects in the time dependence of MF's Mi(t),i = 0, 1,2 and 3, 
whereas the similar functions demonstrate a small-scale fluctuation and a small-amplitude 
oscillation in case of the patient with PSE. 

One can note that the sensors with numbers n = 10, 5, 23, 40 and 53 are specific points 
in the brain core of the patient with PSE. It is interesting to observe a dynamical picture 
for the usual and the nonspecific points at the human cerebral cortex. To this aim, the 
results for the nonspecific sensor with number n = 13 are submitted in Figs. (2) - (6). Figs. 
(2) (for the healthy) and Figs. (5) (for the patient with PSE) present the time dependence 
of the first four dynamical orthogonal variables. One can detect a smooth behavior of the 
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Figure 4: The phase trajectories of the phase points Tij(t) = (Wi(t), Wj(t)), healthy person No. 3 
(sensor number m = 13) for an R/ B combination of the light stimulus. 
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Figure 5: The time dependence of the first four orthogonal variables Wj(i),i = 1,2,. ..3 for the 
MEG signals for the patient with PSE, sensor m = 13, an R/ B combination of the light stimulus. 



variables for Wi(t) in the healthy and, in clear contrast, a sharp or irregular dynamics of 
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Figure 6: The phase portraits created by the phase points ^ij(t) = (Wi(t),Wj(t)) from MEG 
signals for patient with PSE, m = 13, an R/ B combination of the light stimulus. 

Wi(t) in the patient with PSE. Figs. (3) (for the healthy) and Figs. (6) (for the patient 
with PSE) show the construction of the phase space made up by separate phase points. We 
depict stratified phase space for the healthy person and for the patient with PSE. Figs. (4) 
show the nonlinear dynamics of the formation of the phase space by the phase trajectories 
for the healthy. 

Here one can observe the pseudoperiodic orbital movement for the phase trajectory in 
the healthy and the quasi-strange attractors in the patient with PSE. In this context one 
must note the small time scales in the dynamics in the healthy and the larger time scales 
in nonlinear dynamics in the patient with PSE. Figs. (7) depict the time dependence of the 
initial TCF and the first three memory functions Mj(i), i= 1,2 and 3 for the patient with 
PSE. Large scale fluctuations and oscillations are visible in the time dependence of Mj(t), i 
= 0, 1, 2 and 3 in the healthy and small scale deformation are evident in the patient with 
the PSE. 

Figs. 8 show the topographic dependence as a function of SQUID-number of the first 
relaxation parameters Ai for a red-blue (R/B) combinations of the light stimulus for the 
healthy subjects in comparison with the patient with PSE. Similar results one can see for 
a red-green (R/G) combinations of the light stimulus. One can observe a large difference 
of the numerical values of this parameter in the healthy as compared to the patient with 
PSE. The parameter Ai differs on average 6-7 times in the all the sensors. One can note 
a specially strong difference in the data between the healthy and the patient in numerical 
values of parameter Ai especially in the sensors with numbers n = 10, 12, 46, 51 and 56 
(for an R/ B combination of the light stimulus) and n = 10, 12, 46, 51 and 53 (for an R/ 
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Figure 7: The time dependence of the first four MF's Mi(t),i = 0, 1, ...3 for the MEG signals in 
the patient with PSE ( sensor number m = 13) an R/ B combination of the light stimulus. Large 
scale fluctuation of all functions become obvious in comparison with a case for the healthy. 



G combination of light stimulus). Both combination of the light stimulus (an R/B and an 
R/G) yield approximately identical results for the sensor's location. 

Figs. 9 (for an R/ B combination of the light stimulus) and 10 (for an R/ G combination) 
demonstrate the behavior of relaxation parameter Ai for each individual of nine healthy sub- 
jects averaged on all sensor location in cerebral cortex in comparison with the patient with 
PSE. A remarkable difference of Ai appears in healthy subjects, this being approximately 
4-8 times, on average 7 times for an R/ B combination in Fig. 9, and 4, 4 times for an 
R/G combination of the light stimulus in Fig. 10. This difference is a reliable indicator of 
serious destruction in functioning of the human organism with PSE. It is important that the 
behavior of the coefficients of A x specifies the singularities of relaxation mechanisms in the 
MEG's signals. From the physical point of view parameter Ri = |Ai| mimics a relaxation 
rate. We see the drastic distinction in relaxation rates for a healthy person as compared to 
the patient with PSE. This fact indicates the crucial role of the specific relaxation processes 
in the pathological functioning of the human cerebral cortex for the patient suffering from 
PSE. 

One notes that the sensors with the numbers n= 10, 12, 46 and 51 locate specific points in 
the patient's brain with PSE. It is necessary to mention the overall positions of the sensors. 
This indicates that the potential abnormality is not confined to the occipital region but 
distributed over the entire brain. 
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Figure 8: The topographic dependence of the first relaxation parameter Ai for nine healthy subjects 
(upper lines) in comparison with the patient with PSE (lower line) for an R/ B combination of the 
light stimulus. The crucial role of the brain zones with sensor numbers m = 10, 12, 46 , 51 and 53 
is clearly visible. 



V. CONCLUSIONS 

Our discussion in this work make it evident that physiological models typically possess 
a great number of parameters, each with their natural range of variability and uncertainty 
in measurement. Relaxation and dynamic behavior of the system signals can vary widely 
from one chosen set of parameters to another. The presented information-theoretic memory- 
function analysis provides one possibility of extracting interrelations within this complicated 
parameter space. The study of the physical and dynamical boundaries between different 
types of behavior is a necessity, both for a better understanding of brain function and for 
the application of a diagnosing procedure and treatment of patients suffering PSE. 

Control can be applied at preventing the brain from entering an undesirable, pathological 
state such as a seizure [32} . Here we have shown that the parameter space of MEG activity in 
the patient with PSE gives rise to a robust chaotic behavior. In order to study spatiotemporal 
cortical dynamics we need to analyze the global MEG data. In this paper we could cognize 
that the relaxation and dynamic singularities may account for the registration of the relevant 
pathological zones in the human cerebral cortex which are responsible for epilepsy. 

Many natural phenomena can be described by distributions exhibiting a time scale- 



invariant behavior [331. I34J] . The suggested approach allows the stochastic dynamics of neu- 



romagnetic signals in human cortex to be treated in a statistical manner and to search for 
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Figure 9: The mean values of the first relaxation parameter Ai for the whole group of the nine 
healthy subjects (n = 1, 2, 3 ... 9), averaged on the total set of sensors m = 1, 2, 3,... 61 versus 
the patient with PSE (m = 10) for an R/ B combination of the light stimulus. One can note the 
drastic difference of approximately of 4,4 times for the healthy vs the patient with PSE! 



its characteristic statistical identifiers. From the physical point of view, the obtained results 
can be put to use as a possible test to identify the presence or absence of brain anomalies as 
they occur in a patient with PSE. The set of our quantifiers is uniquely associated with the 
emergence of time-scale and relaxation effects in the chaotic behavior of the neuromagnetic 
responses of the human brain core. The registration of the behavior of those indicators dis- 
cussed here is thus of beneficial use in order to detect pathological states of separate areas 
(sensors 10, 5,23, 40 and 53 in our case under consideration) in the human brain of a patient 
with PSE. 

PSE is a type of reflexive epilepsy which originates mostly in visual cortex (both striate 
and extra-striate) but with high possibility towards propagating to other cortical regions 
'3^. Healthy brain may possess an inherent controlling (or defensive) mechanism against 



this propagation of cortical excitations, the breakdown of which makes the brain vulnerable 
to trigger epileptic seizures in patients 36|. However, the exact origin and dynamical na- 
ture of this putative defensive mechanism is not fully known. Earlier we showed [3l[ that 
brain responses against chromatic flickering in healthy subjects represent strong nonlinear 
structures whereas nonlinearity is dramatically reduced to minimal in patients. 

There are other quantifiers of a different nature, such as the Lyapunov's exponent, 
Kolmogorov-Sinai entropy, correlation dimension, etc., which are widely used in nonlin- 
ear dynamics and related applications, see in Ref. In the present context, we find that 
the employed statistical and dynamical measures are not only convenient for analysis, but 
are also suitable for identification of anomalous brain behavior. The search for other quan- 
tifiers, and foremost, the optimization of such measures when applied to complex discrete 
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Figure 10: The mean values of the first relaxation parameter Ai for the group of the nine healthy 
subjects (m=l, 2, ...9) averaged on the whole set of sensors with numbers 1 < n < 61 vs the patient 
with PSE (m= 10) for an R/ G combination of the light stimulus. The distinction between the 
healthy and the patient with PSE amounts up to 8 times! 



time dynamics presents a true challenge. This objective particularly holds true when at- 
tempts are made to identify and quantify an anomalous functioning in living systems. This 
study presents a first stepping stone towards understanding of nonlinear brain processes 
defending against hyper excitation to nickering stimulus by new analysis techniques based 
on non-Markov random processes. 
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